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1 Introduction 

Let T) be the sheaf of analytic differential operators of n variables , . . . , x„. We 
consider a left ideal I oiV generated by {£i, . . . ,£m} which are in the Weyl al- 
gebra D = C{xi, . . . , Xn, di, . . . , dn)- If no confusion arises, we also denote by / 
the left ideal ...,£„} in £>. Assume that Vjlvs, holonomic. It was proved 

by M.Kashiwara that the germs of the fc-th extension group Ext\,[T> j 1 ^ O) form 
a finite dimensional vector space over the field of complex numbers C [Sj ■ We 
note that the vector space is called a fc-th order (cohomological) solution space. 
In [5], an algorithm by which to determine the dimension of the vector space 
was given. In the present paper, we will present an algorithm by which to 
construct a basis of this vector space in a free module over the formal power 
series. In 0, we studied the adapted free resolution oi D/I and an algorithm 
of computing restrictions of D-modules. The algorithm of evaluating the di- 
mension of the germ of the fc-th extension group was an immediate application 
of Cauchy-Kowalevski-Kashiwara's theorem on the restriction of the 2?-module 
T)/! to the origin and the fc-th extension group. In the present paper, we will 
explicitly construct matrix representations of boundary operators of complexes 
appearing in a proof of the CKK Theorem to construct series solutions. 
Let 

be a free resolution of M = D/I. Then, for a left Z?-module N, the vector space 

Ker {UoiJiDiD''^ N) Homc(L>'''+i , iV)) 
Im {llomD{Db-^,N) RomniD''' , N)) 

is denoted by Ext^(M, N) and is called the fc-th extension group. When M and 
N are holonomic Z3-modules, H.Tsai and U.Walther [7| presented an algorithm 
by which to determine a basis of the extension group. Let us consider the 
sheaf of fc-th extension group 6xtf^{'D/I, O), which can be defined in a similar 
way. In this case, O is not holonomic over D, and hence we cannot apply their 
algorithm; we need a different approach. 



When n = 1, it is relatively easy to determine bases of £xtf^{'D/I,0)o, 
fc = 0, 1 where ! = !)■£. Consider the free resolution 

— >V ^ V/I — > 0. 

By applying 7iomi5(-, C')-functor, we have the complex 

0< — d'^d< — 0. 

Hence, we have 

£xt%,{V/I, O) = {/ e C> I £ • / = 0}, £xt\,{V/I., 6) = d/£*d 

following the definition of Ext. Algorithmic methods by which to determine 
bases of the vector spaces above are well-known. Among these, we would like to 
examine in greater detail the method explained in the introductory text book 
by T.Oaku jj. The key in this method is to regard O as an infinite dimensional 
vector space over C and to regard the operator ^ as a linear map. Oaku uses a h- 
function to reduce a problem of infinite dimensional vector spaces into a problem 
of finite dimensional vector spaces. Our method is a natural generalization of 
this method. 

The motivation for the present study is the problem of constructing series 
solutions in Gevrey classes of an .A-hypergeometric system Q]. See also for 
the same problem for Lauricella hypergeometric functions. We hope that our 
algorithm can be applied to this problem. In April 2003, the author discussed 
this problem with F.J.Castro-Jimenez, who presented an explicit computation of 
fc-th order solutions of the hypergeometric system for A = (1, 2). This example 
was exciting and has been the motivation for the present study. 

2 Orders in O 

Let i^T be a field and we consider the ring of the formal power series Oq = 
K\\xi, . . . ,Xn]] in rt- variables. We will omit the subscript of O in the sequel. 
When we apply the results of this section for fc-th order solutions, K is assumed 
to be C. We regard O as an infinite dimensional vector space over K by sorting 
monomials of x by an order. In other words, we identify O with a vector space 
by sorting the coefficients of power series by the order. 

Let us firstly consider when n = 1. If we sort the monomials as 1 
then we have the canonical morphism a from to O as 

i^-:={c=(co,ci,C2,...)|Q eif} ^ K[[x]] 

In this example, the coefficients Ck is sorted as cq, ci, C2, . . . ,. 

When n > 1, we have several natural choices to sort monomials in x de- 
pending on a weight vector and a term order; let w G R" be a weight vector 
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satisfying Wi > and -< a term order in Z"g. We sort monomials in x and 
consequently the coefficients of power series by the order -<w For example, 
when w = (1, . . . , 1) and ^ is lexicographic order, we sort the coefhcients 
as c = (co-.-o, cio-.-o, ■ • • , co-.-oi, C2o - Oj • ■ ■)• this case, we define the canonical 
morphism cr by 

X"°0:={c=(co...o,cio...o,...)|c«ei^} ^ . . . , xj] 

Here, a! = ai! • • • a„!, x°' = x"^ ■ ■ ■ x"" . 

Finally, let us discuss how to encode as an infinite dimensional vector 
space over K. In this case, we use the degree shift vector s G Z'' in addition to 
a weight vector w and a term order -<. Let e^x" be the element of where 
is the z-th standard vector. As in the theory of Grobner basis for Z3 131, we 
define ordu,[s](eia;") — w ■ a -\- Si. To define the canonical morphism cr, we sort 
e^x" lexicographically by ordu,[s], —i and -< for a. The canonical morphism is 
defined as follows 

Here, i = 0, . . . , r — 1, a G Nq and eo, . . . , Gr-i are standard unit vectors. 

We denote by {K"°°Y [s]<„i the image by a^^ of ejx" of which ordu,[s]- 
degree is less than or equal to m. There exist the following natural projections 

T„: iK^°°y[s] (if"-)'' [s]<„ 

Tm'm-- iK^^Y[s]<m' ^ (if"-)'' [s]<„ 

Here, rn' > m and (K"°°Y [s] is the union of {K"°°Y [s]<m for the natural 
numbers m in {K"°°Y . 

We call T the truncation map. 

Example 2.1 Put w = (1, 1) and s = (0,0, —1). We denote 6^x^X2 by k]. 
We consider (if )'^[(0, 0, — 1)]. The elements of which degree is less than or 
equal to 1 are 

[0,0,0], [1,0,0], [2,1,0], [2, 0,1], (degree 0, 4 elements ) 

[0, 1, 0], [0, 0, 1], [1, 1, 0], [1, 0, 1], [2, 2, 0], [2, 1, 1], [2, 0, 2] (degree 1, 7 elements). 
The projection ti^ is a map from if to if*. 



3 Adapted resolution and induced linear map 

We fix a weight vector w in the sequel. We consider a complex of left iJ-modules 
with degree shifts m G Z^, m' G Z«, m" G Z'' 

DP [m] 4 L»9 [m'] 4 D'' [m"] (2) 
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(see |S] on the degree shift). We suppose that the foUowing two conditions are 
satisfied. 

(1) ker(i3) is generated by rows Ai of the matrix A. (2) in(-_^„ „,')[m'](ker(i?)) 
is generated by in(_„, „,)[m'](Ai). When these two conditions are satisfied, the 
complex is cahed adapted at the object I?''[m']. A free resolution is called 
adapted if it is adapted at every object in the complex. The notion is introduced 
in [5] and a free adapted resolution can be constructed by a Grobner basis 
method for a given left D module D/I, a weight vector w, and a term order -<, 
which is used as a tie-breaker. See [H] for an efficient construction algorithm. 

Let us consider the ring of formal power series O = K[[xi, ■ ■ ■ , Xn\], and a 
complex 

DP[s]^ D'^[t] ^ U'[u] (3) 

which is adapted at the middle object D'^[t]. Note that we suppose that D^, 
D^, are sets of row vectors. By applying Homx)(-,C') to the complex, we 
have 

dp[s]^ d'>[t]<^ &[u] (4) 

where O*, (i = p, q, r) are regarded as sets of column vectors. By the canonical 
isomorphism a Q induced by the degree shift, the weight vector and the tie- 
breaking term order ^, we have the following complex of linear vector spaces 



u 



(5) 



Here, and A"^ are block upper triangular matrices with the elements in K . 
The blocks are partitioned by ord^, [-J-degree. The two properties will be key 
ingredients of our algorithm. Since it is block upper triangular, we obtain the 
following complex of finite dimensional vector space by truncating to the degree 
m 

{K^^Y [u]<™ [t]<,„ {K^°°Y [s]<™ (6) 

Let us illustrate the matrix representations of the boundary operators A'" 
and B"^ in terms of boundary operators in (5.2) of 5 , which are used to compute 
restrictions of D-modules. We consider the quotient Ker (A»)/Im (£•) in 
First, we note that 



Ker (A.) 



A.f = 0, f = 



and 









\ 




g = 










\gr/ 


t 
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We denote by the {j, £)-th element of the matrix A. Then the condition 
• f = is equivalent to 

(0) = 0, for all z G N^' and J = 1, . . . , p. (7) 

Define a^g by 

normallyOrdered(9M^,)|_„ = J] ay,a^ a^, G K. 

k 

Here, normallyOrdered(L)|^^g means that 

(1) expand L £ D into normally ordered expression as ^ aakx^d^ 
and then 

(2) replace all Xi by 0. 

In terms of a^^^, the condition ((TJ is equivalent to 

E 4/f ^(0) = 0, for all * e NJJ, J = 1, . . . ,p (8) 

where f^^'^ is the k = (fci, . . . , A:„)-th derivative of fi £ O. Note that Cj^ in JQ) 
is equal to /^"■'(O). Therefore, under the morphism cr, Ker(yl») is nothing but 
the kernel of the matrix defined by {ali)- 

Next, we consider B. Let i?^ be the (j, ^)-th element of B. Define 6^''^ by 
normallyOrdered((3*i3^)|^^g — X^^w^'^- Then the i-th coefficient of the series 
expansion of B^ • gp is expressed as 

Therefore, under the morphism ct, Im is nothing but the image of the 
matrix defined by 

The matrices [alg) and [b^jig) agree with those appearing to compute restric- 
tions of I?-modules (5| Theorem 5.3]. Let us explain this fact. The boundary 
operators A and B induce the following if-linear maps 

where = D/ [xiD + • ■ • + XnD). The cohomology groups of this complex are 
called restrictions. See (Sj for details on computing restrictions. Let us construct 
matrix representations of A and B. We denote by e-' the j-th standard vector. 
The vector space i7(8)£)DP[s] is spanned by {i G Nq, j = 1, . . . ,p), which is 
sent by the hnear map -A to normallyOrdered(9* J2i ^i^^)|x=o — J2k e ^ke^^^^- 
By sorting 9*e^ by the same order to sort the coefficients of power series (we use 
the correspondence d^e^ <-> x'^ei-i), we conclude that the matrix representation 
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of A agrees with the transpose of the A''" in The same assertion holds for 
B. This relation will be used in the proof of Theorem 13. II 
Let M — D/I he a, left holonomic D-module and 

an adapted free resolution associated to a weight vector w. We assume that the 
complex (PJ is a part of this adapted resolution. Let fci be the maximal integral 
root of the 6-function of M associated to the weight {—w,w) 

Theorem 3.1 The truncation map 

KerT,„/(A'^) KerT„,(A'^) 
■ ImT,„,(BT) ^ lu,T^{BT) 

is an isomorphism of vector spaces when m' > m > fci . 

Proof. The matrix t„/(^"^) is block triangular. Define submatrices Afj by 

Submatrices Bf^ are defined analogously. Since JBJ is a complex for any m, 

we have Bf^Af^ = 0, Bf^A^^ + B'^^A^^ = 0, B'^2^'^2 = 0. We note that 
the boundary operators AJ2 and B22 agree with those in the complex (5.3) 
of [S]. Since the complex (5.3) in is exact when m' > m > fci, we have 
ImS^ = Ker A^2- Let {p,q) be in Kerrm'(^"^); we assume that pAji — and 
PA21 + qA'^2 = 0- Then, we can define a natural projection 

i}} : Kerr™/(A^) 3 (p,g) e Ker Afi. 

By utilizing the properties of the matrices A^ , Bf^ stated above, it is easy to 
check that tp induces a well-defined map "0 from KerTm/(A-^)/ImTm'(S-^) to 
Kerrm(^ii)/ImTm(i3ii) that ■0 is injective. It follows from Theorem 5.3 
of El that the dimensions of ^"'"^'"'^^ ^ and ^'""^"'("^^ ) agree. Therefore, we 
conclude that ip — Tm\m is an isomorphism. Q.E.D. 

The next theorem immediately follows from Theorem 13. II 

Theorem 3.2 Retain the notation of the proof of Theorem XH.lX When m' > 
m>ki, fore € Kerrm(A"'"), there exists a vector c' such that (c, c') G Kerrm'(^"^). 
In other words, the linear inhomogeneous equation with unknown vector d 

c • + c' • 74^2 = 
is always solvable for c satisfying c ■ Aji = 0. 
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4 Algorithm 



Put Ai =V/I. We are ready to state our algorithm to construct a basis of the 
d-ih. order solutions. It follows from Theorems 13 . II and 13 . 21 

Algorithm 4.1 Construction of d-th cohomological solution fxtp(A^, O) 
Step 1. Construct an adapted resolution [if^i) of D/I. We assume that the 
resolution is written as © at the degree d. Note that we have assumed A = 

ipd+i, B = ipd- 

Step 2. Let ki be the maximal integral root of the 6 function of M with respect 
to the weight vector {-~w,w). 
Step 3. Obtain a basis of 

Imrfc,(B^) 

as a ii'- vector space. We denote the basis by c^^^ , ■ ■ . , c*^*^) . 

Step 4- By repeating to solve the linear equation in Theorem 13.21 extend the 

vector c*^*) to an infinite dimensional vector in [t] 

Step 5. Output a{c^), . . . , a{cx) G -Rr[[a;i, . . . , as a basis of the solutions. 

Remark. 

(1) In Step 3, our implementation chooses a basis in Kerr^j (v4^)n (imrfe^ (B^)) 
where V denotes the orthogonal complement of the vector space V by the 
standard innerproduct. 

(2) Note that we have to truncate the iteration of Step 4 when we execute this 
algorithm on a computer. 

(3) Let J be a submodule of D^. A basis of £xt^{V^/J, 6) can be constructed 
in an analogous way. 

The proof of the following theorem follows from discussions in Section O 

Theorem 4.2 The set of power series cr(c^''), . . . , cr(c^'') is a basis of 









|f e 6" 


A 













(9.\ 




1' 








\gr/ 





as a K-vector space. In particular, when d — 0, £xt^{M, O) is a basis of the 
classical power series solutions since the denominator is in the expression 
above. 

Example 4.3 Put 

I = D ■ {xdx ~ x{xdx + ydy + 2){xdx + 3), ydy - y{xdx + ydy + 2){ydy + 5)} 
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and consider Ad = D/I. The ideal / annihilates the AppcU function -^2(2, 3, 5, 1, 1, x, y). 
An adapted resolution of M with respect to the weight vector (—1,-1,1,1) is 
as follows. 



D\-l, 0)] ^ D%Q, 0, -1)] Dm\ 



M 







Here, 



where 



xdx — X — X ydxdy — 6x dx — 2>xydy — 6x 
B=\ydy- xy^dxdy - y^d^ - bxydx - Sy^dy - lOy 

631 

ail 0-12 — 1 + xydy + 6a; — 18xy — lllx^y — 135a;ii/^ 

-ydy + 15y + 45a;y xdx — 9a; — 27a;y 1 



631 = -x^ydldy + x^y'^dldy - x'^y'^d^d^ + xy^d^d^ + ^x'^ydl - Ix'^ydxdy + Qxy'^d^dy - Sxy'^d^ 

15x^ydl+6x'^y'^dxdy-9xy^d^+A5x^yd^+45x^y'^dxdy-27x'^y^dxdy-27xy'^d^ + 

270x^ydx - 135x'^y'^dx + Ubx'^y'^dy - 216xy^dy + 270a;\ - 270xy'^ 

and 

Oil = xy^dxdy - xy'^dl + y'^dl + bxyd^ - 7xydy + B>y'^dy - 5y+33a;t/^9y+60a;t/+ 

162a;2y25^ + l^hxy'^dy + 315a;=^y - 27Qxy'^ - 2565x'^y'^ - 2025xy^ - 5265x^y'^ - 
6075a;2y3, 

ai2 = -x^dl - Sxydy + 3x - ISx'^ydx - 9x'^ydy - Mx"^ + 27xy - lUx^ydx - 
ISbx'^y'^dx - 27x^y'^dy - 27x^y + 1053x^y + 1701a;\2 + SlbQx^y"^ + 3645a;\3. 
Since the 6-function is s, the number fci is equal to 0. The dimension of 

£xt'^{M,d) is 1 and 



E 



(3)m(5)ri 



min 



!n! 



is a basis. 

The dimension of £xt]y{M, O) is equal to 2 at the origin. Let us construct 
a basis of this vector space. The induced linear maps on the space of truncated 
power series ri(A^), 



ri(A^) = 



ri(B^ 



[B"^) are 
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Submatrices standing for the underlined elements are equal to to{A^), to{B^). 
A basis of the vector space ^^^"I^^t-I is (1,0, 0, —5), (0, 1, 3, 0). Here the mono- 
mials BiX^y'^, which is encoded as k], are sorted as [[0, 0, 0], [1, 0, 0], [2, 1, 0], [2, 0, 1]]. 
Then the 0-th approximation of a basis of series solutions is 

(l,0,-5y),(0,l,3x) 

This solution can be extended to the 1-th approximation by solving the linear 
inhomogeneous equation in Thcorem l3.2l for m ~ and m' = 1 as 

(1 + lOy, 0, -by - Abyx + 60?;V2!) 

(0, 1 + 6a;, 3£ - ZQx^ /2\ + 27yx) 
The set of indices is sorted as 

[[0, 0, 0], [1, 0, 0], [2, 1, 0], [2, 0, 1], [0, 1, 0], [0, 0, 1], [1, 1, 0], [1, 0, 1], [2, 2, 0], [2, 1, 1], [2, 0, 2]] 

Repeating the procedure for T^i^)^ T2{B'^), we obtain the following basis 
of second approximate solutions 

1601775000 2 574824600 n 

( l^-10y^ H yx+180y^, 

^ ^ 10611803 10611803 

574824600 278581950 ^ 
10611803 ^ 10611803 ^ ' 

2 38669400 . 39525750 . ^^ 

-5y - 45yx + eOy^ yx^ y^x + 7020y'^) 

^ ^ 10611803-^ 10611803-^ ^ ' 

280878030 2 344820564 ^ 2 344820564 2064283056 2 

( x"^ H yx, l + 6x-|-72x^H yx H y^, 

^ 10611803 10611803 10611803 10611803 

Q 2^07 o-iVR 3 , 15126426 2 , 13920984 ^ ^ 

3x - 36x + 27yx - 2376x H yx ' 

10611803-^ 



10611803 



Here, we set x'y^ = x'^y^ / (ilji). 
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